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Abstract. This paper is a study of the effects of anisotropic matter sources 
on the qualitative evolution of spatially homogenous cosmologies of Bianchi 
type VIII. The analysis is based on a dynamical system approach and makes 
use of an anisotropic matter family developed by Calogero and Heinzle which 
generalises perfect fluids and provides a measure of deviation from isotropy. 
Thereby the role of perfect fluid solutions is put into a broader context. 
*~ j ■ The results of this paper concern the past and future asymptotic dynam- 

ics of locally rotationally symmetric solutions of type VIII with anisotropic 
matter. It is shown that solutions whose matter source is sufficiently close to 
being isotropic exhibit the same qualitative dynamics as perfect fluid solutions. 
C*J , However a high degree of anisotropy of the matter model can cause dynamics 

to differ significantly from the vacuum and perfect fluid case. 
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1. Introduction and motivation 



For spatially homogenous (SH) spacetimes the Einstein-matter equations for a 
large variety of matter sources reduce to an autonomous system of ordinary dif- 
ferential equations in time. Thus the mathematical theory of dynamical systems 
can be applied to gain insights into the qualitative behaviour of SH solutions. This 
{y-^ ■ approach has been used in mathematical cosmology, e.g. to address questions rele- 

vant for observational cosmology, in particular concerning the role the Friedmann 
solutions play in the more general context of SH cosmologies that are not spatially 
isotropic in general. On the other hand the interest in SH models is nourished by 
the believe that the dynamics of SH cosmologies towards the initial singularity is 
crucial for the understanding of the behaviour of more general spacetimes close to 
singularities; cf. [1] and references therein. 

Less is known about SH solutions with matter sources more general than perfect 
"k>(( i fluids. Calogero and Heinzle have developed a matter family naturally generalis- 

ing perfect fluids that contains large classes of anisotropic matter sources and is 
suited for a dynamical system analysis. It includes a measure of the deviation from 
isotropy and thus allows to investigate the role of perfect fluid solutions in the more 
general context of solutions with anisotropic matter. The SH cosmologies consid- 
ered were of Bianchi type I, and locally rotationally symmetric (LRS) types I, II 
and IX; cf. [2, 3]. While dynamical system analyses with specific anisotropic matter 
sources have been carried out before, the approach by Calogero and Heinzle can 
be regarded as a first step towards a systematic study of the effects of anisotropic 
matter to the qualitative dynamics of SH cosmologies. 

This paper is concerned with the analysis of SH cosmologies of LRS Bianchi 
type VIII with anisotropic matter. In section 2 the basic features of the anisotropic 
matter family are stated, and the state space and evolution equations for LRS 
type VIII are given. The dynamical system analysis is performed in section 3, where 
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some technical details on the analysis of the flow at infinity are contained in the 
appendix. The results are given and discussed in section 4, where the main result is 
formulated in theorem 1 and corollaries 1 and 2. As a small byproduct, the results 
cover the future asymptotics for perfect fluids with p — wp and w € (— |, 0), which 
might fill a little gap in the literature. Section 5 is concerned with an extension of 
the formalism to treat Vlasov matter dynamics with massive particles. 

Part of the material needed in sections 2, 3 and 5 has already been presented 
in [3] or [4] to analyse LRS Bianchi types I, II and IX. At these points in the text, 
only the crucial steps and results are quoted from there. 

2. The LRS Bianchi type VIII setup 

In a frame (d£, a) 1 , w 2 ,w 3 ) adapted to the symmetries, an LRS Bianchi class A 
metric has the form 

4 g = -di® dt + gn(t) uj 1 ® uj 1 +g 2 2(<) (w 2 <8> a) 2 + a) 3 <E)uj 3 ). 

In type VIII, da) 1 = —^€jkin i: 'uj k A w', with [n* J ] = diag(— 1, 1, 1). Greek indices 
denote spacetime components while Latin indices label spatial components w.r.t. 
the adapted frame. The metric will be subject to the Einstein equations — without 
cosmological constant — in geometrised units (8irG = c = 1), cf. [3, Eq 2]. 

2.1. The anisotropic matter family. For a perfect fluid that is non-tilted with 
respect to dt, the components of the stress-energy tensor are [T M „] = diag(/?, p,p,p), 
where p and p denote the energy density and pressure of the fluid. It is an isotropic 
matter model since the eigenvalues of \T i j\ are all equal. When p — wp with 
w = const, the fluid is said to obey a linear equation of state. 1 

The matter models considered in this paper form a family of models generalising 
perfect fluids with linear equation of state. This family of models is described in 
detail in [3, section 3]. In the following, a brief description tailored to the present 
purposes is given: 

The components of the stress-energy tensor are [T M „] — &ia.g(p,p\,p2,P2) where 
the energy density p > and the isotropic pressure p :— YliPi/^ obey a linear 
equation of state p = wp, with w = const. Defining the dimensionless rescaled 
principal pressures as Wi :— Pi/p it is clear that w — JV Wi/3, which in turn 
implies that w\ and W2 are not independent once w is given. As a consequence 
of the Einstein equations and some basic assumptions on the matter family given 

22 ( 1 . 

in [3, section 3], w-i (and thus w\) is a function of the quantity s := y? u 6 (0, ^J. 
Clearly, s gives a measure of anisotropy of the spatial metric components while W2 
encodes the anisotropy of the matter content." The simple example of isotropic 
matter corresponds to w\ — W2 — w. Also, note that s — > and s — > ^ correspond 
to 522 — > oo (or gn — > 0) and gn — > oo (or 522 — > 0) respectively, i.e. to a singular 
metric. A basic assumption is that the limit of w\ for s — > ^ coincides with the 
limit of u>2 for s — > 0, i.e. that wi (5) = 102(0). Hence one can define an anisotropy 
parameter 

(1) P:=2 W - W2{0) 

1 — w 

that provides a measure of deviation from an isotropic matter state in the extremal 
cases where the metric is singular. Note that fi — corresponds to matter models 
that behave like a perfect fluid close to singularities. 

The analysis in section 3 will show that the qualitative dynamics of LRS Bianchi 
type VIII solutions does not depend on the whole function wi but merely on the 



For example, w = corresponds to dust and to = i to radiation. 
In [3] Wi(s) is called anisotropy function and denoted by u(s). 
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value 1^2(0) where the metric is singular. Thus any two matter models of the 
anisotropic matter family that share the same parameters w and /3 also share the 
same qualitative dynamics, even if they have different functions Wi(s). Accordingly, 
P serves as the parameter to investigate the influence of the anisotropy of the matter 
on the dynamics, even though it gives a precise measure of the anisotropy only close 
to singularities. Therefore, in the context of the qualitative analysis of the dynamics 
of solution that follows, a specific class of matter models is simply characterised by 
a pair (u>, /3) in the parameter space 

P:= (-|,l) x R. 

Here w is restricted to (—h, l) since the primary interest is in matter models that 
obey the standard energy conditions [5, p 218-220]: The weak energy condition 
corresponds to m; > —1. The strong energy condition requires the weak energy 
condition to hold and w > — |. The dominant energy condition is \wi\ < 1. 
Therefore, by (1), for the energy conditions to hold, one has to restrict /3 to 
[max(— 2, — j3jj), l] ; cf. [3, table 2 in section 3.4]. The dominant energy condition 
is only marginally satisfied when /3(w) takes the boundary values. The parameter 
space P and the subset for which the energy conditions are satisfied is illustrated 
in figure 3 together with the bifurcation lines which will be explained in section 3. 

2.2. The Einstein equations as a dynamical system. LRS Bianchi types VIII 
and IX share the same evolution equations which have been derived in detail in [3] . 
This subsection gives a brief outline: 

Due to spatial homogeneity, the Einstein equations for LRS type VIII with the 
above matter source are a constrained autonomous system of ordinary differen- 
tial equations in t for the components of the spatial metric g^ and the extrin- 
sic curvature fey. This system is then written in terms of quantities that are 
standard cosmological parameters and/or bring the equations into suitable shape. 
These quantities are the Hubble scalar H := — h(k l i + 2fc 2 2), the shear variable 
(j + := ^(fc 1 ! — fc 2 2) and the quantity mi :— J (711/522- Finally these are divided by 
the dominant variable D := -J ' H 2 — 1/(3.922) > to obtain the normalised variables 
(fr D) E +l Mi) := {H,a + , mi )/D. 3 

The resulting representation of the Einstein equations for anisotropic matter 
filled LRS Bianchi type VIII spacetimes is then given by the dynamical system 

-(l-H 2 D )(q-H D X+) 



(2) 



Mi 



-(2 - q)H D V+ - (1 - H 2 D )(l - Z 2 + ) + ^ + 3 ft ( W2 ( s ) _ w ) 
Mi (qH D - 4E+ + (I - #£)£+) 



which represents the evolution equations, and the constraints 

M 2 I 
(3) Q + Ei + -i = l and 1 -#£, = -—- < 0. 



The former is the Hamiltonian constraint, while the latter follows directly from 
the definition of D. ft := ^y > denotes the normalised energy density and 
q := 2S 2 _ + 1+ 2 3 "' £1 the deceleration parameter. The prime denotes derivatives 
with respect to rescaled time; (•)' := jjjj( - )- Note that the definition of Mi and 
the second constraint imply that s can be regarded as a function of Hp and M\, 

s = (2 — 3 M i D ) . Therefore (2) is a closed system once u>2(s) is prescribed 
through the 'equation of state' of the anisotropic matter. 



In [3] this normalisation has been taken in lieu of the Hubble normalisation [6, section 5.2] 
mainly because it yields a compact LRS type IX state space. Although this is not the case for LRS 
type VIII, the dominant normalisation still has favourable properties. For example the type III 
form of flat spacetime is a fixed point solution in this formulation; cf. table 1 . 
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Figure 1. The state space Ayin and its boundary subsets. 

2.3. The state space AVm. The LRS type VIII state space is determined by 
the constraints (3): Since J7 > and Mi > by definition, it follows from the 

Hamiltonian constraint that £+ < 1 and £+ + -jj- < 1 . The inequality 1 — H D < 
imposed by the second constraint implies that the state space is the union of two 
disjoint sets: Hd > 1 corresponds to positive H and hence to forever expanding 
universes, while Hd < — 1 corresponds to forever contracting universes. However, 
since (2) is invariant under the reflection (t,Hr>,H+) — >• — (i,i?£>,S+) it suffices to 
restrict to the expanding case. Accordingly the LRS Bianchi type VIII state space 
is defined as 

~H D -\ 

H D G (1, oo), £+ G (-1, 1), M 1 G (0, ^12(1 - Z%)) 



X 



VIII 



M x 



which forms the tunnel-like structure depicted in figure 1. The boundary subsets 
of AViii arc Vyin, 6m and 5jj, which correspond to tt = 0, Mi = and Hd = 1, 
respectively. LRS type VIII vacuum solutions thus lie in Vyin while solutions in iSj 
and Bin are of Bianchi types II and III respectively; cf. [3, sections 9.2 and the first 
remark in 10.1]. S^, which corresponds to Hd —* °o, is of course not a boundary, 
but one can think of it as boundary in a compactified version of the state space; cf. 
appendix A. 

In this formulation LRS type IX cosmologies are subject to the same evolution 
equations (2). However there is a sign change in the second constraint of (3) leading 
to a different state space X\x for which Hd G (— 1, 1). Figuratively speaking Xix 
builds the tunnel connecting Ayiii with the second disjoint LRS type VIII set. 
Hence both Ayin and Xix share the boundary 5j. Cf. [3, section 9]. 

There are two challenges in connection with the analysis of the dynamical sys- 
tem (2) in Xym. First, Ayin is not compact. Hence one has to perform a careful 
analysis of the 'flow at infinity'. In the present case it will be shown that there do 
not exist orbits that emanate from or escape to infinity. Second, the dynamical sys- 
tem (2) does not extend to the line Ci in dX^ui for that Hd = 1 and M\ = since 
s(Hd, Mi) has no limit when C\ is approached from AVin- 4 The analogous problem 
occurs in the context of LRS type IX, where it was overcome by introducing another 
set of coordinates that give regular access to this part of the boundary by 'blowing 
up' Cf, cf. [3, section 10.2]. This method can be adapted to the LRS type VIII case 
by applying the same coordinate transformation. However the corresponding state 
space is again different: 

2.4. The state space 3Viii- The coordinate transformation used to analyse Ci is 
(4) 1-H D = 2rcos8, M\ = 3rsin6>, S + unchanged, 



However s(Hp, Mi) has limits when C\ is approached from B\n or S»; cf. sections 3.2 and 3.3. 
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Figure 2. The state space 3-Viii and its boundary subsets. 



where r > and 9 £ [f ,tt]- Thereby, Ayiii is transformed to the state space 



JViii 



r £ 0, 



4(1 



i6 



,£+e(-i,i) 



which is depicted in figure 2. From (4) one has the following correspondences 
between subsets of AViii an( i 3Vm : 



sfux!, 



(5) Ayiii ~ Jfym, Vvm ~ ^ viii i ^ni ~ B lu UXj, <Sj ~ Sj UAfj, £j ~ A"i. 

iSoo corresponds to a line at infinity in the context of DVin- Furthermore, (4) defines 
a diffcomorphism between A'viii\>Ci and 3Vrii\^i> from which it follows that the 
flows in these sets are topologically equivalent. In particular, 

(6) Xyjji S yvni, Vvm = V^jjj, Bm = Bf u and S$ = S^ . 

In contrast, the coordinate transformation performs a blowup of the line C\ to the 
two-dimensional set X\ defined by setting r = 0. It can be identified with the LRS 
Bianchi type I state space; cf. [3, section 10.2]. With (4), s can be regarded as 
a function of 9 alone. Therefore s has a limit as r — > 0, which implies that the 
dynamical system, when expressed in the coordinates (r, 6, £+), extends regularly 
to Xi. 

Finally, note that since s(9) is a bijection on [f,7r] one can as well choose s as 
coordinate instead of 9; hence the two labels on the axis in figure 2. 

3. The dynamical system analysis 

Whenever possible the analysis is carried out in the coordinates (Hd,T, + , Mi), 
which is in all sets except X\. However the final results presented in section 4 have 
to be interpreted in the context of the state space 3-Viii since the system (2) is not 
regular in AV111, which prevents a complete global analysis in the original state 
space. 

Since the matter parameters w and (3 enter the evolution equations everywhere 
except in the vacuum boundary and at infinity, cf. sections 3.1 and 3.5, their 
values determine the qualitative properties of the flow. For instance there are 
fixed points that only occur in the state space iff (w,j3) is in a certain subset of 
P. Similarly, fixed points may have different local stability properties depending 
on (w, 0). The curves /3(w) £ P dividing P into these subsets corresponding to 
qualitatively different dynamics shall be called bifurcation lines. These bifurcation 
lines will be found in the subsequent subsections and are plotted in figure 3. 

3.1. Analysis in Vyiii- The dynamical system in Vvm is obtained by setting 
fi = in (2) and using (3): 



(7) 






(1 - Hl){H D - 2S + )S + 
(1 - S5.) (2 + (1 - YR ) + (H D - E H 
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Figure 3. The bifurcation diagram in the parameter space P. The 
shaded region (including P D ) marks the subset for which the energy 
conditions are satisfied. 



There are three fixed points in Vvni, T := [1, -1] T , Q := [f , 1] T and D := [2, 1] T . 
The eigenvectors and eigenvalues of the linearisation of (7) at the fixed points 
determine their local stability properties. One finds 



T 



"6" 
12 




'1 0" 
1 


, Q: 


' 2 " 

-4 




1 0" 
1 



and D 



1 



so T is a source, Q is a saddle repelling in i/p direction and D is a sink in Vvni- 

It is proven in appendix A.l that there are exactly two more fixed points at 

infinity, T^ = [oo, — 1] T and Qoo = [oo, 1] T . The stability properties then follow 

directly from (7): First, E' + > in Vyin- Second, for Hd > 2 (and thus at infinity), 

H' D = for E + = 0. Hence X^ and Q^ play the role of saddles for the flow in 

Vvni given in figure 4. 

— y . . — 

To interpret this as flow in V VIII , note that the fixed point T in Vvni corresponds 

— y 
to the closure of the line Tj, — > T$ in V V in- Since Tj, and Tp turn out to act as a 

— y 

source and a saddle in V vm in all cases respectively, cf. figures 5 to 17, the orbits 
in Vym have the form T\, — > D. 



3.2. Analysis in Bui. The dynamical system in Bui is obtained from (2) by setting 
Mi = (•«• s = 0) and using (1): 



(8) 






-(1-fTl,) (2- |(1 -«;)(! 



ij£)S H 



-(1 - Y? + ) ((1 - F|,) + |(1 - W )(ff B S + + /?)) 



-.:f(H D ,'E + ) 



Here and henceforth the notation follows the pattern P : 



lyx, V2] where A^ and Vi denote 



the i eigenvalue and eigenvector of the linearisation of the dynamical system at P. 
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H D 



FIGURE 4. The flow in Win- 

This system has at least three and at most five fixed points in Bui depending on 
(w, j3) G P, namely 

2+3f3(l-w) 



T, 



r i 




it 




r?" 




\ 1 1 


-i 


,Q\,-= 


i 


,D = 


i 


i R\> '■— 


.-0. 



and P 



^/(l-3w) 2 +6fS(l^w) 

l+3w 

^/(l-3w) 2 +6/3(l-w) 



Clearly R, is a fixed point in Bm (and different from Tj,, Q\,) iff j3 G (—1,1). The 
conditions H D \ P > 1 and £+|p < 1 entail that P is in Bm iff /3 > j3\,(w) := ^- 
and (w,/3) $■ F D , where P D refers to the small D-shaped subset of P bounded by 

cf. figure 3. Under these conditions the square root in 



. . -l±^-3+(l-3w)2 



the coordinates of P is automatically real and S+|p > 0. 

The eigenvectors and eigenvalues of the linearisation D/(iJo,E_| 
and R\, determine their local stability properties. One finds 



at Tb,Q\,,D 



XL 



Qb 



6 




'1 0" 


3(l-u-)(l-/9)_ 




1 


2 




'1 0" 


3(1-uj)(1+/3)_ 




1 



D : 



R ? 



3/3(1 



-3 

w) 



6w 



l-3tu 



3/3 2 (l— w) + l+3iv+2/3 



l-2w+P(l-w) 
1 

3^ 2 (l-tu)+5+3tu+4^ n 
1 1 



Therefore, in Bm, T[, is a source for /3 < 1 and a saddle repelling in -ffp direction 
for /? > 1, Q\, is a saddle repelling in iJc direction for /3 < —1 and a source for 
/3 > — 1, -D is a sink for /3 < /3\,(w) and a saddle attracting in i?o direction for 
/3 > /3\,(w) and i?i, is a sink \/(w,(3) G P^ and a saddle attracting in S + direction 
VK-1</3<1)GP\P D . 

The eigenvalues X± of D/|p are given in terms of trace and determinant by 

(9) 
where 



(10) 
(11) 



detD/|p 



2A± = trD/|p ± v/(trD/|p) 2 -4detD/|p, 
_ 9(1 - w)(3^ 2 (l - w) + 1 + 3w + 2(3) (/3(1 - w) - 2w) 



trD/|p = — 



(l-3w) 2 +6^(l 
3(l-tw)(l + 2/9) 



u>) 



a/(1 - 3w) 2 + 6/3(1 - w) 

which are real whenever P exists in B\u because the square root in (11) appears in 
the coordinates of P as well. Further, trD/|p < and detD/|p > whenever P 
is in Bm. 1 ' Hence the real part of (9) is always negative, whether the eigenvalues 
are real or complex, so P is a sink in Bm whenever present. The eigenvalues are 
real for /3 G (A (u>), /??(«;)]; cf. figure 3. 



To see the first inequality, note from figure 3 that > 



when P is in Bin. To see the 



second inequality, note that setting the middle and last factor in the numerator of (10) to zero 
corresponds to the bifurcation line /3±(w) and /3\,(w) t respectively. 
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In order to perform the analysis of the flow in the full state space <%vin m 
section 3.5, it is also necessary to know the local stability of P in the direction 
of Mi. From (2) one has (lnJV/i)'|p = — 3£+|p, which is negative so that P is 
attracting in the direction of M\. This implies that P is a local sink in AViii- 

It is shown in appendix A.l that there are exactly two more fixed points at 
infinity, Too = [oo, — 1] and Qoo = [oo, 1] . Their stability properties then fol- 
low straightforwardly from (8): For ifp, sufficiently large, f(Hjj,'E+) « [— J?| ) E+, 
fff,(l - X^)] T . Hence, £^ > and H' D = for £+ = 0, which implies that T^ 

and Qoo play the role of saddles for the flow in Bm. 

For the cases where P does not exist in 23m, the information suffices to draw 
the corresponding qualitative flow diagrams. When P exists in Bm, periodic orbits 
encircling this fixed point could in principle be present. However numeric investi- 
gations strongly suggest that this is not the case. In any case, the main results of 
this paper are not affected by this open question; cf. section 4. 

The resulting qualitative dynamics in Bm in dependence of (w, /3) £ P is depicted 

— y — y — 

in figures 5 to 17 as flows in B m together with the flows in Sa and Xj. The latter 

two subsets will be discussed next. 

3.3. Analysis in Sp Sj is the common boundary of the state spaces Ayin and Xix 

of LRS types VIII and IX. It has been analysed in detail in [3, section 10.1]. The 

bifurcation lines f3 = ±2, f3$(w) and (3Uw) result from this analysis. The results 

—y , , —y — 

are depicted in section 4 as flows in Sa together with the flows in B lu and X\. 

3.4. Analysis in X\. The flow in X\ has been analysed in detail in [3, section 10.2]. 
The only new bifurcation line in figure 3 resulting from this analysis is (3 = which 
corresponds to the local stability properties of the Friedmann point F. However 
the qualitative dynamics in X\ is also associated with the bifurcation lines [3 = ±1 

— y 

and p — ±2. The results are depicted in section 4 together with the flows in B m 
and 5jj . 

3.5. Analysis in AVm and Soo. The full system (2) does not have any fixed 
points in Aym- Furthermore, as shown in appendix A. 2 there are only the already 
know fixed points Too = [oo,— 1,0] T and Qoo = [oo,l,0] T at infinity. Their sta- 
bility properties then follow straightforwardly from (2): For Hjj sufficiently large, 
rhs(2) w [-#!,£+, #1,(1 - £+), -23 r £,E + Mi] T . Hence £' + >0,H' D ==0 for £+ = 
and M[ = for S + = 0, which implies that T^ and Qoo play the role of saddles for 

the flow in AViii- Viewing infinity as the boundary S^ in a compactified version 
of the state space, T^ is a source and Qoo is a sink in <Soo, hence all orbits in S^ 
are of the form T^ — > Qoo. 

An important consequence is that there is no orbit in Xym that emanates from 
or escapes to infinity. Hence each orbit in Ayin lies in a compact subset. The 
following lemma concludes this section. It will also be used to localise the a- and 
cj-limit sets [G, p 99 Def 4.12] in the next section; cf. [6, p 91 Def 4.7] for the term 
'invariant set': 

Lemma 1. Let 7 be an orbit in 3Vm- Then both, a("f) and coi'j), is a non-empty, 
compact and connected invariant subset of D^VTII- Furthermore, a(j) C Su U X\ 
and w(7) C B m U X\. 

Proof. First, since each orbit in Ayin lies in a compact subset, so does by (4) each 
orbit in 3Viii- The first statement of the lemma then follows from [6, p 99 Thm 4.9]. 
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Second, consider the function Z5 : Ayin U Vviii — > K given by 

Z 5 := 5- > with Zr = — — Z 5 < 0. 

(#f, - 1)1 2H D 

One can check that Z$ |s + =n=o = and Z'^'\s + —q—o < in Ayin, which means 
that Z§ is strictly monotonically decreasing along the flow in Afyin- Hence the 
monotonicity principle [6, p 103 Thm 4.12] implies that the limit sets lie in SihUiSj. 
Moreover, since Z$ — > = inf(Zs) when one approaches B\\\ and Z5 — >• 00 = 
sup(-Zs) when one approaches iSj, the monotonicity principle implies 01(7) C 5j 
and u (7) C Bm for all orbits 7 in Ayiii- In the context of 3Vin this means that 
a(7) C Sjf U #1 and a; (7) C 1^ U #1; cf. (5). □ 

4. Results and discussion 

Finally all information obtained in section 3 can be collected to identify the a- 
and w-limit sets of generic orbits 7 £ 3Vm, and hence the past and future asymp- 
totic dynamics of generic LRS Bianchi type VIII cosmologies, for all qualitatively 
different anisotropic matter cases: 

4.1. Identification of the limit sets. Figure 3 shows that there are thirteen 
qualitatively different cases.' These cases are listed in figures 5 to 17, which show, 

for each case, a plot of the bifurcation diagram where the corresponding subset of 

— y — — y 
P is shaded and a representative flow diagram in the set S* U X\ U 8 m , on which 

the limit sets lie by lemma 1. Solid (dashed) lines represent generic (non-generic) 

orbits in the respective boundary subset. A filled (empty) circle on a fixed point 

indicates that it attracts (repels) orbits in the respective orthogonal direction. s 

Heteroclinic cycles and networks are represented by thick lines. The axes with two 

labels emphasise the diffeomorphisms (G) and the bijection between s and 9; cf. 

section 2.4. 

By lemma 1 , the limit sets are non-empty, compact and connected invariant sub- 
— y 
sets of |V viii- This leaves fixed points, heteroclinic cycles and heteroclinic networks 

as candidates for the limit sets in figures 5 to 17. 

A single fixed point is an (a-) w-limit set of generic orbits 7 € JVni iff it (repels) 
attracts orbits from a three dimensional neighbourhood in 3A/m; c -g- a (source) 
sink. These can easily be identified from the figures. In the cases of figures 9 to 1G, 
fixed points are the only possible generic limit sets. Furthermore, these cases exhibit 
just one source and sink respectively, which implies that these are the past/future 
attractors [G, p 100 Def 4.13]. 

Figures 5 to 8 and 17 show a heteroclinic cycle or network each, and hence 
additional candidates for limit sets. In figures 7 and 8 the cycle is the past attractor 
since there is no other candidate for a generic a-limit set in 5u U<li, and as proven 

in section 3.5 0(7) is non-empty. Note that since the cycle is not in B ul U X\ it 
cannot be an w-limit set. The same line of arguments shows that in figures 5 and 6 
the past attractor is a subset of the heteroclinic network. 

In figure 17, the heteroclinic network could be an a-limit set in addition to Tj,. 
Analogously, in figures 5 and 6, the heteroclinic cycle dX\ could be an w-limit set 
in addition to D and P, respectively. The nontrivial task of directly tackling the 
stability properties of these structures is omitted here. They are however taken 



The bifurcation lines j3T(w) and K(w) only specify if the fixed points Cj and P are local 
stable nodes or foci in Sj and B^ u respectively, which is irrelevant for identifying the limit sets. 

The colour coding of Too and Qoo makes only sense in the context of Aym- The fixed points 
in dX\ are not colour coded since the orthogonal directions are represented in the figures. 
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Fixed points corresponding solution 



ype 


vacuum 


I 


/ 


I 


• 


I 




I 




II 

TTT 




in 


/ 



Tj,,Tjj Taub Kasner 

Q\,,Q$ non-flat LRS Kasner 

R\,,R$ a type I anisotropic matter solution 

F flat Friedmann 

Cj generalisation of Collins-Stewart 

P generalisation of Collins (VI_ i) 

D type III form of flat spacetime 

Table 1. The exact solutions corresponding to the fixed points. 



into account as further candidates for generic limit sets in these three cases. Note 
that figure 1 7 is the only one of these three cases where the energy conditions can 
be satisfied; cf. section 2.1 and figure 3. 

4.2. The main results. Finally one arrives at the main theorem and its corollaries: 

Theorem 1. The a- and co-limit sets of generic orbits 7 £ JVin of the dynamical 
system (2) with parameters (w,/3) £ P, is given as stated in the captions of figures 5 
to 17. These describe the past and future asymptotic dynamics of LRS Bianchi 
type VIII cosmologies with anisotropic matter of the family defined in section 2.1. 

Proof. Cf. section 4.1. □ 

The exact solutions corresponding to the fixed points are summarised in table 1 . 
These have been given in [3, appendix A] for all fixed points but D. For D, one 
just needs to insert the coordinates into [3, Eq 99]. The corresponding isotropic 
cases to these solutions can be found in [6, section 9.1]. 

Vacuum and perfect fluid cases. LRS Bianchi type VIII vacuum solutions are 
past asymptotic to T\, and future asymptotic to D. 

Generic LRS Bianchi type VIII solutions with a non-tilted perfect fluid where 
p = wp and w £ f-|,l) are past asymptotic to Jj,, future asymptotic to P for 
w £ (—|,0) and to D for w £ [0, 1). 

The vacuum part follows from section 3.1; cf. figure 4. For the perfect fluid 
cases, recall from section 2.1 that they correspond to the line j3 — in P. They 
are therefore contained in theorem 1 as special cases. 9 The only bifurcation line 
intersecting (3 = is (i\,{w) at w — which yields the two qualitatively different 
perfect fluid cases, see figures 12 and 13, respectively. The statement on the future 
asymptotics for w £ (— -|,0) might fill a little gap in the literature. The other 
results are known: The past asymptotics has been given in [7, table 4 and figure 4] , 
the future asymptotics for vacuum in [8, Prop 8.1] and the future asymptotics for 
perfect fluids with w £ [0, 1) in [9, Thm 3.1]. 

Comparison with the cases of theorem 1 shows that the dynamics with anisotropic 
matter differs from the vacuum and perfect fluid cases in figures 5 to 8, 14 and per- 
haps 17. The strongest implications of this represent the main results of this paper 
stated in the following corollaries: 

Corollary 1 (past asymptotics). The past asymptotic dynamics of generic LRS 
Bianchi type VIII cosmologies with anisotropic matter can differ significantly from 
that of the vacuum and perfect fluid cases. In particular, the approach to the initial 
singularity is oscillatory in the cases depicted in figures 5 to 8 and perhaps 1 7. 



The fact that j3 = is itself a bifurcation line, which is related to the stability of the Friedmann 
point, is not relevant for the asymptotic dynamics. 
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Corollary 2 (future asymptotics). The future asymptotic dynamics of generic 
LRS Bianchi type VIII cosmologies with anisotropic matter can differ significantly 
from that of the vacuum and perfect fluid cases. Note in particular the cases of 
figures 14 and perhaps 5 and 6. 

For completeness it is necessary to elaborate on the remark at the end of sec- 
tion 3.2 concerning the possible occurrence of periodic orbits around P in Bf u : 
Such orbits are not observed numerically. However, their occurrence could change 
certain details in theorem 1, adding further cj-limit candidates in the cases where 
P is in Bf u . On the other hand, corollary 1 would be completely untouched and 
corollary 2 even strengthened by the occurrence of such periodic orbits. 



4.3. Physical interpretation of the results. An interesting observation is that 
there is a neighbourhood of the line f3 = in P for which the corresponding asymp- 
totic dynamics is identical to that of the perfect fluid cases. Hence the perfect fluid 
solutions are robust under small perturbations of a vanishing anisotropy parameter. 
For the past asymptotics to differ from that, f3 even needs to be so high that the 
dominant energy condition is only marginally satisfied. 

The specialisation of the results to the cases for which the energy conditions are 
satisfied is obtained by restricting (3 to [max(— 2, — y^jj), l] ; cf. section 2.1. This 
corresponds to the shaded region in figure 3, and thus rules out the cases of figures 5 
and 6. All other cases contain values (w, (3) for which the energy conditions are 
satisfied. 

However, it should be pointed out that the cases represented by figures 7 and 8 
only satisfy the energy conditions for /3 = 1, and the case of figure 17 only for 
w > h and (3 = —2, i.e. for (w, (3) in one-dimensional subsets of P for which the 
dominant energy condition is only marginally satisfied. Consequently, in these cases 
the energy flow is necessarily lightlike close to the singularity. A special example is 
given by collisionless (Vlasov) matter with massless particles, which as shown in [ ', 
section 12.1] falls into the class of matter models with (w,f3) = (i, l). 

In any case, the statements of corollaries 1 and 2 remain true under the restric- 
tion to matter satisfying the energy conditions. 

This concludes the main results on the anisotropic matter analysis. Section 5 
is concerned with an extension of the present formalism to treat LRS type VIII 
dynamics with Vlasov matter with massive particles, which does not a priori fit 
into the anisotropic matter family considered so far. 
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Figure 5. 0(7) C network(X|,,Tjj,(2jj,(2i,)- 0^(7) = D, possibly dX\. 
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Figure 7. 0(7) = cycle(T b ,T|j,Qtt,Q b )- w (t) = -P- 
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Figure 8. 0(7) = cycle(T b , T h Q h Q b ). w( 7 ) = D. 
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Figure 9. 0(7) = T b . w( 7 ) = £>. 
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Figure 10. 0(7) = T b . u(y) = P. 
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Figure 11. 0(7) = T b . w( 7 ) = P. 




2 
1 




r /3 
/ 




/ 




















1 


1 



4 


r 





2 


,s 








12 


Mi 





~ 





7T 


1 



/ir 



Figure 12. 0(7) = T b . 01(7) = D. 
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Figure 13. 0(7) = T b . w( 7 ) = P. 
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Figure 15. 0(7) = T b . 01(7) = L>. 
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Figure 16. 0(7) = T b . 01(7) = L>. 
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Figure 17. 0(7) = T b , possibly network^jQ^F). 1^(7) = D. 
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5. Extension of the formalism to treat Vlasov matter dynamics with 

massive particles 

It was stated in section 4.3 that Vlasov matter with massless particles falls into 
the class of matter models described in section 2.1 with parameters (w, /3) = (4, 1). 
This is shown in [3, section 12.1]. However, as stated there, for Vlasov matter 
with massive particles the relation p — wp is non-linear. In [4], Calogero and 
Heinzle extended their formalism to treat Vlasov matter dynamics with massive 
particles in the case of LRS type IX. Since types VIII and IX are analogous in this 
formulation, and in particular share the same evolution equations, this extension 
of the formalism can be adopted to LRS type VIII without difficulties: 

Following [4, section 3], the dynamical system representing the LRS type VIII 
Einstein- Vlasov system is given by the system (2) and the additional equation 

(12) I' = 2H D l(l - I). 

The additional variable I := \_,°, 9 ' , 1/3 G (0, 1) corresponds to a length scale of the 
spatial metric. The rescaled principal pressures are functions of I and s. Their full 
expressions are given in [4, Eq 13] from which it follows that 

(13) (tfl,j8)|l=0 = (|,l) and K/?)| i=1 = (0,0). 

The LRS type VIII state space for Vlasov matter dynamics is given by Xym x (0, 1). 

The dynamical system can be analysed as follows: By (12), I is strictly mono- 
tonically increasing along orbits in AViii x (0, 1). The monotonicity principle thus 
implies that 0(7) C AV111 x {0} and oj(j) C Afyin x {l}- 10 Hence, the search for 
the limit sets can be restricted to these boundary subsets: 

From (13), the flow in AV111 x {0} is equivalent to the flow for anisotropic 
matter with parameters (ijl); cf. figure 8. Hence the cycle (T^Xj, Q$, Q\,) is the 
past attractor of the flow in AV111 x (0, 1), since this cycle is the only candidate 
for an a- limit set for generic orbits. Also from (13), the flow in AV111 x {1} is 
equivalent to the flow for anisotropic matter with parameters (0,0), i.e. to that for 
dust; cf. figure 12. Hence D is the future attractor of the flow in Ayin x (0,1), 
since this fixed point is the only candidate for an w-limit set for generic orbits. 

I = corresponds to Vlasov dynamics with massless particles; cf. [4, section 3]. 
Hence the result states that Vlasov matter with massive particles behaves like 
Vlasov matter with massless particles asymptotically to the past, and like dust 
asymptotically to the future. 

Appendix A. Fixed points at infinity 

To find the fixed points of (2), (7) and (8) at infinity, a method presented in [10, 
3.10] is adopted in a slightly modified way: 






A.l. In two dimensions. Consider a two dimensional system 

where P and Q are polynomials of degree n and m in Hp respectively, satisfying 
n < m + 1 . One can write this system as 

(14) Q(H D ,Z + )dH D -P(H D ,Z + )dZ+=Q, 

where however the information about the direction of the flow is lost. Next the 
Hd coordinate is formally compactificd by projecting the state space onto the 



One can use the same arguments as in section 3.5 and appendix A for each I 
hypersurface to prove that all orbits can be trapped in a compact subset of R 4 . 
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Figure 18. Projection onto the 'Poincare cylinder'. 



'Poincare cylinder' as illustrated in figure 18. The corresponding coordinate trans- 
formation (iT D ,£ + ) -> (X,E+,Z) : X 2 + Z 2 = 1 is given by H D = |-. Applying 
this to (14) and multiplying by Z m+2 yields 

(15) Z m+1 Q(X/Z, £+)dA - Z rn+2 P{X/Z 1 S+)dS+ - XZ m Q{X/Z, S+)dZ = 0, 

which defines a flow on the 'Poincare cylinder'; cf. [10, p 267]. Points at infinity in 
the original state space correspond to points on the 'equator' X = 1, Z = on the 
'Poincare cylinder'. Furthermore, evaluating (15) with X = 1, Z = gives the flow 
on the 'equator', which corresponds to the flow at infinity, 

(16) (Z m Q(X/Z, £+)) \ x=ltZ=0 dZ = 0. 

Note that the first two terms of (15) vanish since they are at least proportional to Z. 
This is not true however for the third term of (15). From (16), for {Z m Q)\x=i,z=a ^ 
it follows that dZ = 0, which corresponds to trajectories through a regular point 
on the 'equator', where the sign of (Z m Q)\x=i,z=o determines the flow direction. 
Fixed points on the 'equator' correspond to solutions of (Z m Q)\x=i,z=o = 0; cf. 
[10, p 268 Thm 1]. Solving this equation for Q in the context of (7) and (8) 
yields S + = ±1 for the fixed points at infinity in Vyin and B\u respectively, i.e. 
Too := [oo, -1] T and Qoo := [oo, 1] T . 

A. 2. In three dimensions. The generalisation to three or more dimensions is 
straightforward; cf. [10, p 277 ff]. Consider a three dimensional system 



~P{H D) T. +) M{) 
Q{H D ,V + ,M X ) 

_R(H D ,Z+,M 1 ) 

where P, Q and R are polynomials of degree I, n and m in Hd respectively, satis- 
fying n + l>l<m+l. One can write this system as 

Q(H D , £+, M X )&H D - P{H D , £+, Mr)dS + = 

(17) 



£+ 

Mi 



R{H D ^+,M X )&H D - P(H D ,'£ + ,M l )dM l = 0, 

where however the information about the direction of the flow is lost. Again the 
Hjj coordinate is formally compactified by means of a projection on the 'Poincare 
cylinder', i.e. by a coordinate transformation (iJu,E + ,Mi) — > (X, Y, + , M X ,Z) : 
X 2 + Z 2 = 1 given by Hd = ^r- Applying this to (17), multiplying by Z m+2 
and Z n+2 respectively and evaluating the resulting expressions at the 'equator' 
X = 1,Z = yields 

(Z m Q{X/Z, £+, M 1 ))\ x=lz=Q dZ = 

(Z n R(X/Z, £+, M 1 ))\ x=lz=Q dZ = 0. 

Trajectories through a regular point on the 'equator' corresponds to dZ = 0, where 
the flow direction is determined by the signs of (Z m Q(X/Z, S + , Mi)) \ x=1 z=0 
and ( y Z n R(X/Z, S + , Mi)) I _ . Fixed points on the 'equator' correspond to 
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solutions of the system of equations {{Z m Q)\x=i,z=o = 0, {Z n R)\x=i,z=o = 0}. 
Solving this for Q and R in the context of (2) yields £+ = ±1, Mi = for the fixed 
points at infinity in Ayin, i.e. Too := [oo, — 1,0] T and Qoo := [oo, 1,0] T . 
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